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Abstract 

Strongly coupled J\f = 4 supersymmetric Yang-Mills plasma at finite temperature 
and chemical potential for an R-symmetry charge undergoes a second order phase 
if^ . transition. We demonstrate that this phase transition is of the mean field theory type. 



We explicitly show that the model is in the dynamical universality class of 'model B' 
according to the classification of Hohenberg and Halperine, with dynamical critical 
exponent z = 4. We study bulk viscosity in the mass deformed version of this theory 
in the vicinity of the phase transition. We point out that all available models of bulk 
viscosity at continuous phase transition are in conflict with our explicit holographic 
computations. 
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1 Introduction 

According to gauge theory/string theory correspondence of Maldacena [1] maximally 
supersymmetric A/" = 4 SU{N) Yang-Mills (SYM) theory is dual to string theory on 
AdS^ X S^. In the planar limit {qym — ^ 0, A^ — )> oo with A = Qym-^ ^^P^ fixed) and 
for large t' Hooft coupling A ^ 1 the strongly coupled SYM is described by classical 
type IIB supergravity on AdS^^ x S*^, making it essentially soluble. The value of this 
holographic duality is that it can provide explicit tests of various phenomenological 
models invented to describe the dynamics of strongly coupled systems. The focus 
of this paper is the application of gauge/gravity duality to the transport properties 
of strongly coupled gauge theory plasma in the vicinity of the second order phase 
transitions. 

In [2] it was argued that the only model for the critical behavior of the bulk viscosity 
in strongly coupled systems at continuous phase transitions not in confiict with explicit 
holographic computations was that proposed by Onuki [3]. Specifically, Onuki's model 
predicts that close to the phase transition the bulk viscosity scales as 

Coc |t|-^'^+", (1.1) 

where 

t = ^-l, (1-2) 

is the reduced temperature, v and a are the usual static critical exponents of the 
continuous phase transition, and z is a dynamical critical exponent. In this paper we 
would like to definitely answer the question as to whether or not (1.1) is realized in a 
strongly coupled gauge theory plasma with a holographic dual. 

Our starting point is the best studied example of gauge theory/string theory duality, 
namely that of A/" = 4 SYM plasma. This theory has an 5*0(6) ~ S'f/(4) R-symmetry; 
thus one can turn on three independent chemical potentials (one for each of the f/(l)'s 
in the Cartan subalgebra of the R-symmetry group). It is well known that A/" = 4 SYM 
plasma at finite temperature T and for a single U{1) R-symmetry^ chemical potential 
/i undergoes a second order phase transition^ [4, 5, 6]. Moreover, recently [7], the 



^This is not the diagonal U{1) of the SU{A) R-symmetry. 

^ As we show below, some of the static critical exponent first computed in [5] and since then widely 
used in the literature are incorrect. This issue could be traced back to the fact that the hyperscaling 
relation between static critical exponents is violated in this theory. 



conductivity aq of this gauge theory plasma was shown to be finite on the critical line 

critical V ^ 

As a result, the authors of [7] argued that the dynamical universality class of A/" = 4 
SYM plasma is that of 'model B' according to classification of Hohenberg and Halperine 
[8], with the dynamical critical exponent 

z = A-r], (1.4) 

with 7] being the anomalous static critical exponent. In this paper we confirm the 
identification made in [7], and compute z for the A/" = 4 SYM plasma. 

Unfortunately, we can not use the A/" = 4 SYM plasma directly to test Onuki's 
prediction for the scaling of the bulk viscosity in the vicinity of the phase transition 
(1.1) — conformal invariance of the theory guarantees that the bulk viscosity must 
vanish for arbitrary chemical potential and the temperature. Thus, we need to deform 
the theory in such a way that we break the scale invariance. The simplest deformation 
one can consider is to give mass M to fermions of A/" = 4 SYM. If 

M < T,r^tical, (1.5) 

it is sufficient to work to order O ( ^ ) . Although not necessary, one can think about 
above deformation (to the order specified) as that corresponding to deforming A/" = 4 
plasma to A/" = 2* plasma [9, 10, 11]. 

The paper is organized as follows. We being section 2 with presenting the effective 
action for the holographic dual of A/" = 4 SYM plasma, deformed by a dimension-3 op- 
erator. We study this gravitational model at finite temperature and chemical potential 
to order O ( ^ ) : we (numerically) determine the background geometry, discuss the 
holographic renormalization of the theory, and compute the equilibrium thermodynam- 
ics. We present a highly nontrivial test on our analysis by demonstrating that the basic 
thermodynamic relations are satisfied. Appendix A contains necessary technical de- 
tails. In section 3 we study hydrodynamic fiuctuations in charged gauge theory plasma 
and their holographic dual — the lowest quasinormal modes of the (mass-deformed) 
Reissner-Nordstrom (RN) asymptotically AdS^ black holes. We derive the speed of 
sound and the sound wave attenuation coefficient in a generic charged plasma. We 
find that the attenuation coefficient is sensitive to both the shear and the bulk viscosi- 
ties of the plasma, as well as the plasma conductivity. Interestingly, the dependence 



on a conductivity in the attenuation coefficient arises only at the fourth order in the 
parameter breaking the scale invariance, i.e., ^. We derive (coupled) equations for 
the quasinormal modes in the deformed RN black hole background and explain how to 
decode from their spectrum the speed of sound waves and the bulk viscosity of the dual 
plasma. We point out the computational difficulty in the numerical analysis, intrinsic 
to finding the sound channel quasinormal modes in charged black hole backgrounds, 
and present a new method of computing the corresponding quasinormal modes. Ap- 
pendix B is used to explain this new method in a simple setting of A/" = 4 SYM plasma 
at finite temperature, but zero chemical potentials. As a highly nontrivial test of our 
analysis, we show that the speed of sound obtained from the thermodynamic analysis is 
in excellent agreement with the one extracted from the dispersion relation of the lowest 
sound channel quasinormal mode. In section 4 we discuss dynamical critical phenom- 
ena in A/" = 4 SYM plasma. We review how the dynamical susceptibility can be used 
to compute static anomalous critical exponent rj and the critical exponent u associated 
with the divergence of the correlation length in the vicinity of the transition, as well as 
the dynamical critical exponent z. We explain how the problem of finding the sound 
channel quasinormal mode can be adjusted to extract the dynamical susceptibility. We 
emphasize why the latter analysis can not be performed in the hydrodynamic limit. 
We present results for the critical exponents for A/" = 4 SYM plasma and demonstrate 
that these exponents are robust against the mass deformation of A/" = 4 plasma. We 
summarize and interpret all the results in concluding section 5. 

2 Holographic dual of mass deformed Af = 4 plasma at equi- 
librium 

Effective five-dimensional action describing the holographic dual to A/" = 4 SYM de- 
formed by an operator Oa of dimension A takes form^ 

^5 =-^ f d'^y^ C 



167r( 

' (2.1) 



^The supergravity gauge coupling is chosen so that the asymptotic AdS:^ radius is one. 



where 5C is a mass deformation 

6C = -\{dxf-'^x' + 0{x'). (2.2) 

As usual, the mass of x is related to the dimension A of the corresponding operator in 
the dual description 

A(A - 4) = m^ (2.3) 

In what follows we focus on A = 3 {rn^ = —3) deformation. 

The non-normalizable component A of x near the (asymptotic) AdS^ boundary is 
related to the coupling M of operator O^ deforming the Lagrangian Ccft of the A/" = 4 
conformal fixed point: 

Ccft -^ Ccft - MO3 , M oc A , (2.4) 

where the precise definition of A and the relation between M and A (up to an irrelevant 
c-number normalization of O3) will be established later. Notice that the identification 
of (2.4) with the dual holographic action (2.1) can be established only to order A^, and 
thus x^ in (2.2). The latter fact is emphasized by an O(x^) term in (2.2). The reason 
why this is so is best illustrated with the M = 2* example of the holographic gauge 
theory/string theory correspondence [12]. The duality studied in [12] is one of the few 
examples where it is possible to match exactly the gravitational parameters with the 
corresponding gauge-theoretical ones. In this case, a simple deformation of the type 
(2.4) leads to a complicated potential for a supergravity scalar x — nonetheless, up 
to order M^, the dual supergravity deformation is unambiguously fixed by the scaling 
dimension of Oa, as in (2.2). 

We are interested in the critical phenomena in A/" = 4 plasma (and its massive 
deformation) at finite temperature and chemical potential — thus as far as M ^ T, in 
particular for T = T^., if is sufficient to work with the effective holographic description 
(2.1), (2.2). 

2.1 Background 

In this section we set up our notations for describing the background geometry of 
A/" = 4 SYM plasma at finite temperature and U{1)r chemical potential, deformed by 
O3 operator to quadratic order in its coupling. 
Consider the following ansatz 

dsl = -cf dt^ + cl dx^ + C3 dr"^ , A^ = AS^, (2.5) 



where Cj = Ci{r) , A = A{r), = 0(r), x = x(^)- We find it convenient to introduce a 
new radial coordinate x as follows 



1 — X = — 

C2 



(2.6) 



so that a; — )> 0+ corresponds to an AdS boundary and x — )■ 1_ corresponds to a regular 
Schwarzschild horizon. Further, we introduce 

1/2 



^% C2 = g 



H 



(2.7) 



with g = g{x), H = H{x). We would like to construct background geometry pertur- 
batively in x, in other words we parametrize the background as 

^ = ^o + A2^2, H = Ho + X^H2, A = Ao + \'A2, X = Axi, (2.8) 
where the coefficient of the non-normalizable mode of x, ^-c, A, is introduced so that 

;^i = x^/^ + 0{x^^^) , as x-^0+. (2.9) 



To order (9(A°) we find 



go = 13: 



Ao 



(3VT 



K 



K 



m 1 



K 



H. 



2(1 + k) + {2x - x^)k^ + K^x{2 - x){2 + xk){2{1 + k) - xk) 



(2.10) 



2{1 + k) 

where the two constants {(3, k,} are related to the temperature T and the chemical 
potential /x of the R-charged black brane 



2ttT 



/3 



K + 2 



A=0 



Note that the ratio 



^k{1 + k) 


2nT 




/i 


A=0 



/i 



/S 



A=0 



VT+K 



v^ + 



K 



attains a minimum at /t(A = 0) = 2, at which the black brane undergoes 
order phase transition [4, 5, 6]. 
■ To order 0{\^) we have 



(2.11) 

(2.12) 
a second 



=Xi + 



„ H^{1 + k) + H^il + k)+ HI(2k^ -k + 2k^x^ - 1 - Axk^) -1-k 



{HI - 1){HI + l){x - l){l + n) 
H^K'^m^{Hl{K^x^ - 2xK^ + HI-2k + k'^ -2) + K+l + H^k) 
(1 + KfiHl - l)\Hl - 1)2 



X'l 



Xi 



(2.13) 



■ To order (9(A^) we have 

=g'i + Cn 92 + C12 H'^ + C13 A^ + Cu (x'l)' + C15 ^2 + Cie H^ + C17 x? , 
=H'^ + C21 (7^ + C22 H'^ + C23 A'2 + C24 (x'l)' + C25 ^2 + C26 ^2 + C27 X? , 
=Ai + C31 (7^ + C32 H'^ + C33 A'2 + C34 (72 + C35 H2 , 

where the coefficients Cij are collected in Appendix A. 



(2.14) 



2.2 Holographic renormalization and the boundary stress-energy tensor 

In this section we carefully perform the holographic renormalization of (2.1) and extract 
the thermodynamic quantities, subject to the following constraints: 

■ the background geometry of the gauge theory dual to the gravitational action (2.1) 
is i?3,i; 

■ since (2.1) is expected to be valid only to order (9(x^), we perform holographic 
renormalization only to this order in x, but exact in 0. 

Let r be the position of the boundary, and S*^ be the Euclidean gravitational action 
on the cut-off space 

\\mS^E = SE, (2.15) 

where Se is the Euclidean version of (2.1). Explicitly, using equations of motion. 



SI 



levrG 



dr 



5 Jr^ 



d ^^JHeCe 



1 



SvrG. 



dM 

4 



levrG 



dr 



5 Jri^ 



d^iVhC 



dM 



dX 



dMs 



dr 



rh 



2 1-'' 
C1C2C2 

C3 



SttG. 



dMa 



d'a X 



,2^/ 1 



C1C2C2 
C3 



(2.16) 



where rh is a position of the regular Schwarzschild horizon of (2.5), He is the induced 

metric on the boundary, and 

1 N^ 

Notice that for a regular Schwarzschild horizon, the horizon contribution in (2.16) 
vanishes. Besides the standard Gibbons-Hawking term 



5, 



GH 






1 [ci4]' 



d'^, 



(2.18) 



dMs 



we supplement the combined regularized action {S^ + Sgh) by the appropriate bound- 
ary counterterms which are needed to get the finite action. These boundary countert- 
erms must be constructed from the local metric and {(f), x} scalar invariants^ 



S 



counter 



1 



^nG, 



5 JdMr, 



C^ W ^S ttl + "2 X +"3 (111' 



o{x' 



(2.19) 



ai are constant coefficients of the counterterms which are determined by the require- 
ment of having a finite action. The counterterm oc 0:3 is finite, but is required to 
insure the first law of thermodynamics [15] and the supersymmetry (at extremality) 
[14] of the single-charge AdS^ RN black hole. Thus, the finiteness of the renormalized 
Euclidean action 



Ie = lim I Se + Sgh + S' 



counter 






\Ie\ < 00, 



constraints 



3 1 1 

«i = 7: , "2=0' cts 



2 ' 8 6 

These values of ai agree with those reported in [17, 18, 14]. 

The boundary stress energy tensor T^'^ is obtained from the variation 
action 

C C^ I C I cf^ounter 

Otot — ^E + ^GH + 'J , 



(2.20) 
(2.21) 

of the full 

(2.22) 



with respect to the (Minkowski) boundary metric rj^i, = c "^ h 



flU 



rp^U Q 






' 87rG, 



,2 



-e^"" + Qh^"" - 2<^ ai + a2 x' + "3 (ln0)' \ h^"" + 0{x 



where 



e^'^ = - {v^n'^ + v^n^) , e = Tr e^^ 



(2.23) 



(2.24) 



2.3 Asymptotics of the background geometry 

Within the background ansatz (2.5)-(2.7), it is straightforward to derive from (2.1) 
equations of motion for A,H,g,x and solve these equations perturbatively near the 



^In principle, there could be a (finite) contribution from the gauge field A^. Such contribution was 
shown in [14] to vanish. 



boundary x — )■ 0+, and near the horizon y = 1 — x —)■ 0+. These equations enjoy exact 
scahng symmetry 

{A, g, H, x) ~ (A A, A g, H, x) , (2.25) 

for any constant A. 

2.3.1 The boundary x — )■ 0+ asymptotics 
We find 

A =ao + ai x^/2 + f^ _ 3/,^ _ LA ai x + O (x^/^ i^ 3,) ^ (2.26) 



191/9 /a?/iiCn 7cn 1 /iiCn a? , 



+ 0(x3/2ln2x)') , 



24 ^ V 96^^ 1728 4 24 8c/^^ 

(2.27) 

if =1 + /ii x^l'' + f^ + ^ Vi ^ + ^ (^'^' In ^) , (2-28) 

X = Co x^/^ (1 + ( cio + |j Inx J x^/2 + C (xlnx) J . (2.29) 

The general solution is determined by 6 parameters 

ao, ai, ^1, fl'o, Co, Cio j , (2.30) 

of which one is the temperature (it can be thought roughly as a scaling parameter in 
(2.25)), and another other two {ao, cq} are the coefficients of the non-normalizable 
modes related to the f/(l)i? chemical potential fi = a^ and the mass-deformation scale 
M oc Co respectively. 

The remaining parameters {ai, hi, cio} correspond to the charge density p (conju- 
gate to the chemical potential), and the expectation values of the dimension-2 {O2) oc 
hi and the dimension-3 {O3) oc coCio operators. To establish precise correspondence, 
we need to introduce a 'universal RG scale' — a universal (in a sense of being indepen- 
dent of any scales in the gauge theory: the temperature, the chemical potential, and 
the mass-deformation scale) radial coordinate: 



r = C2 



'^^^"'0Hw4)^"^^^(^"-')- ^'■''' 



The precise values (up to a c-nuinber normalization) of various mass scales and oper- 
ator expectation values (VEV's) can then be obtained as the coefficients of the non- 
normalizable/normalizable modes of the corresponding gravitational fields with a radial 
dependence given by (2.31): 



2 / \ 2 



M = cq — -^-—p, , p = --^-pr «i — T^TJ ) (C'2) = hi — -^ 

3i/2/,y2 SttGs \3^/^hl^y \3^/^hl^\ 

(2.32) 

The factor of {—StiG^)^^ in the definition of p is inserted to insure consistency of the 
thermodynamics of the dual gauge theory. Notice that the expectation value {O3) is 
determined only up to order 0{cq) — this is related to the fact that the effective action 
(2.1) is defined only to this order as well. 

2.3.2 The horizon y = 1 — a; — )■ 0+ asymptotics 

We find 

A = 4y' + 0{y') , (2.33) 

^ = 4 + 92 y' + 0{y') , (2.34) 

^ _,H (4(/^o^)^(4)^ + {h^,)\a^,nc^,Y + Kh^,f{a^,? + ^9^^^^? - 1) , 

3(^o")^(/^S(co")^ + 4) ^ 

+ 0{y') , 



(2.35) 



The most general solution is determined by 5 parameters 

U, 9l 9^ hi A . (2.37) 

Note that once the temperature, the chemical potential and the deformation scale 
are fixed, we have 5 + 6 — 3 = 8 parameters, which is the correct number to uniquely 
determine the solution. 



10 



2.4 Thermodynamics of the background geometry 

Given (2.32), the thermodynamics of the mass-deformed theory must be studied under 

the constraint 

31/2/^1/2 

Co = - -^M, (2.38) 

9o 

for a fixed M. 

All the thermodynamic quantities are given to order 0{cq) or (9((co)^). 

The temperature T and the chemical potential /x of the background geometry is 
given by 

mcH + ^g^' 

The entropy density s is 



/9^^n2 l"-0l'-0J ^^)Kyo) ,, ^ /^ qQ\ 



^_Aor_ ci _ 1 27r{g^)'{h^,f/' 



4G., 46-, 



j/->0+ 



^TtGs UK)' - 1)3/2 • 



(2.40) 



The regularized Euclidean action (2.20) (up to the space-time volume factor) has to 
be identified with the Gibbs free energy density Q 

n = ([ d^A~' i, = ^(-4, + 'Mh^lf^, (2.41) 



\JdM, J SttCsV 9/1? 288/1? 

From (2.23), the energy density e and the pressure P are given by 



P 



8nG5 \3hl 288/1? 



8ttG5 \9hl 288/1? 



(2.42) 



Notice that Q = —P, as required by the Minkowski space-time thermodynamics. Fi- 
nally, from (2.32), the charge density is given by 

87rG5 3hi 
While not automatically satisfied, we expect the basic thermodynamic relations 

Q = e — Ts — fxp , de = T ds + fi dp , dP = s dT + p dp , (2.44) 

to hold. 
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2.5 Thermodynamics to order (9(A°) 

The background geometry to this order is given by (2.10). Thus, we find (we do not 
need cio) 

ao, a,, hi, go, coj = (-7=- -2'/'/?, ^4f=' /^' o) ' (2-45) 



and 



4, ^0, ^2, K, c'o] = { ^/Y+^^^ + 2y ^' °' ^^ + ''^'^'' ° ' ' ^^-^^^ 



leading to (we use (2.17)) 



P 



27r(l + kYk^^T^N^ 



{K + 2y 





(« 


: + 2)4 


27iT 


- v^ + 


2 



(2.47) 



It is straightforward to verify the basic thermodynamic relations (2.44). 

From (2.47) we see that - achieves a minimum at k = k^ = 2, corresponding to 
the critical temperature T^ = \/2fi/TT and the critical chemical potential jic = 7iT/^/2. 
Introducing 

t = — -l, u = l--^ =^ n = , 2.48 



we find 



27ArV^ A , 8 , _ 16v/2 „,/2 , 6^ 



„ _ > 21N^-TV /. 4 lev's ,„ 14 , _,^ , ,„, , 



(2.49) 



(2.50) 



K =K±(t) = 2± 4V2 t^/2 + 8 t ± 5v/2 t^/2 + 4 t^ + C (t^/^) 

=/t±(/i) = 2 ± 4^2 ;U^/2 + 8 /i ± 7^2 ^^/2 + 12 ^^ + C (/i^/^) , 

where the signs in both expressions correlate. Thus for a given temperature t there are 
two thermodynamic phases of the system, with r2_ being the stable one. Specific heat 
of the Q+ phase is negative, and thus this phase is thermodynamically unstable. In [16] 
it was argued that in the absence of the chemical potentials in holographic examples of 
gauge theory/string theory correspondence the thermodynamic instabilities show up 

12 



as instabilities in the propagation of sound waves. In RN plasma of interest here there 
are no instabilities in the sound channel; rather, we find that in the Q^ phase the two- 
point correlation function of the order parameter (the charge density) in the vicinity 
of the critical point oscillates, instead of exponentially decaying. Additionally, the 
dynamical relaxation time in the f2_|_ phase is negative, which explicitly demonstrates 
the instability of this phase (see section 4). 
For i7_ phase we find: 



-Kl^) 



oc 



P 



dn^ 



dfi 



oc 






oc 



+ri/2 , 



dp 



OC 



J, djj, 



oc +p 



i/^oc+t^/^ 



(2.51) 



oc +/i 



-1/2 



OC +t 



-1/2 



rp dp? 

in the vicinity of the phase transition. 

Recall that in the classic theory of static critical phenomena one introduces critical 
exponents 

(a,/3,7,<5,//,^) . (2.52) 

Once we identify the charge density p of the RN plasma with the order parameter, and 
the chemical potential p with the external (control) parameter we can read from (2.51) 
the following 4 critical exponents 

1 1 1 



(a,/3,7,(5) 



2 '2'2 



(2.53) 



Under a single scale hypothesis in a continuous critical phenomena only two static 
critical exponents are independent — there are 4 scaling relations: 



a + 2/3 + 7 = 2, 7 = /3(5 - 1) = z/(2 - r]) , 2 - a = vp , 



(2.54) 



where p is the number of spatial dimensions of a critical system. Notice that given 
(2.53), the first two scaling relations in (2.54) are satisfied. One can use the rest of the 
scaling relations to determine the remaining static critical exponents: 



V 



relation 



V 



(2.55) 



relation 



As we will see in section 4, the values of {i',i]) as in (2.55) are incorrect. We point 
only that the value of r] in (2.55) already raises suspicion: we consider here the critical 
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phenomena in RN plasma in the strict 't Hooft hmit, i.e., as A^ — )■ oo — in this hmit 
we expect a mean-field criticality, leading to 



V 



= 0. (2.56) 

mean— field 



2.6 Thermodynamics to order (9(A^) 

Since we are going to be interested in the deformation of the critical phenomena in 
(2.49) at order C(A^) (or more precisely C(M^), see (2.38)) we have to solve (2.13) 
and (2.14) in the vicinity oi k = 2. Thus, we further expand 

OO OO CO 

g, = P J](«: - 2)" gi-^ , H, = J](«: - 2)" H^^ , A^ = /3 ^(/. - 2)" 4"^ , 

n=0 n=0 n=0 

OO 

Xi = $^(/.-2)"xS"^ 

n=0 

(2.57) 

Clearly, k = 2 is not a singular point of (2.13) and (2.14), thus (2.57) leads to a series 
of smooth ODE's for {g2 , -f^2 ? 4 5 Xi }■ These ODE's must be solved with the 
following boundary conditions: 

a;-^0+: ^f ^ , if ("^ ^ , A^^ ^ , xt^ -^ ^^o ^''\ 

X — 7> 1_ : (72 ~^ const , if 2 ~^ const , A2' — > , Xi ~^ const . 

(2.58) 

The normalization of the non-normalizable modes of A2 and Xi near the boundary 
implies from (2.39) and (2.32) that 

It is straightforward to construct asymptotic solutions for {(72 > -^^2 ? ^21 Xi }? 
subject to the boundary conditions (2.58) — with obvious modifications they take the 
form of (2.26)-(2.29) and (2.33)-(2.36). Thus we find: 
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for UV parameters (2.30): 



ao 



ai = -(3V2 + \^I3 ^{k - 2)" aj" 



* n=0 



/^i = T4^ + A'E('^-2r/ii"), ^o = /5, co = A, (2.60) 

3vl + /t 

oo 

Cio = ^(fi;-2) 



n=0 

,(n) 
n=0 



^10 ' 



for IR parameters (2.37): 

oo 

+ \'P Y.^K - 2)" a5(") , ^o' = /3 + A^/3 E('^ ' 2)" ^S 



„ oo oo 



(/t + 2)v/r+K 



n=0 n=0 



^,^ = A^/3 5^(«: - 2)" gf-^ , h', = (1 + «:)^/^ + A^ J^^^ - 2)" /^J^") , 

?i=0 ji=0 

oo 

n=0 



where A is given by (2.59). 

Notice that at each order (n) there are eight parameters 

(n) z,(") (") h{n) h{n) h{n) vh(n) h{n) 



(2.61) 



r,(^> /i*-"^ ^i-"^ z,"'-"^ ^"i-"^ ^"i-"^ /i"^"^ ^"t-"^ (0 ao\ 

which are uniquely determined by solving ODE's for {5'2 ? H2 , A2 , Xi }^ subject 
to the boundary conditions (2.58). 

For the first two orders — n = 0, 1 — we collected these coefficients in Table 1. 

Given (2.59), (2.60) and (2.61) we can compute following (2.42) and (2.43) thermo- 
dynamic potentials as a series in {k — 2) and to order (9 ( ^ j . Recall that the basic 
thermodynamic relations (2.44) are not automatically satisfied. Rather, they impose 
algebraic constraints on (2.62). We find: 
■ at order n = 0: 

= 7v^ 3'/' hf^ - ^VS hf^ - \V2 t'^ (cT^r + -V2 gf^ - lOV^ g^^'^ 



+ sVQ aT + af , 



(2.63) 
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n 





1 


aS") 


-0.1655555689119159 


-0.013078142827910066 


hP 


0.11928166889409134 


0.018530943268730934 


An) 
'-10 


-0.23071409088829617 


-0.01240068268619082 


h{n) 
^2 


-0.03250808756595229 


0.005503802347302805 


h{n) 

do 


0.021412846345792092 


0.0015030136731377015 


h(n) 
92 


- 0.007908684866787086 


-0.0009615292197374708 


uHn) 


0.1485225060823788 


0.042145056457606175 


h{n) 
^0 


0.7464562054847809 


-0.0129139123278239 



Table 1: Coefficients of the normalizable modes of the background geometry. 



= -12^3 hf - -V2 3^/3 (cJW)2 + 14^2 3^/3 /^J(o) + ^^ + 32v^ g^ 



h{0) 



24 



18v^ g'/^ 



W^cf^+at 



,{o) 



24v^ a5(°) 



(2.64) 



= - Av/2 cSS + 1^2 cSJ + 4^3 hf^ - ^73 /.« - ^-^Vg af 
16 16 4 3 



18 



Iv/2 3^/^ /.;« + ?v^ 4^°^ - 6v/2 4^^^ - ?v/2 .^^ ' 



6 



9 



^V2 3V3(cr^)^ + a«, 
(2.65) 



at order n 



= -i af^ + |v^ /.f ) - ^72 gr + ^^^2 3^3 (e;(°))^ 
25 
"6 



+ f^^,S(n?v/2,f) 



8v/6 4^^^ - 10^2 



% 



h{i) 



36 ° 

7v^ 3^/=^ hi^'^ - ^4V3 h^P 



-iv^3i/='c;(°)cS(^) + a« 



(2.66) 



= 14v^ 3^/3 hT^ - -V2 3^/3 



^(0)^(1) 



11 



80 



2 >^ ^ ■ 48 

_ 8v/2 g,"(°) + -V2 oT + - V2 c^'il + 32v^ -^"(^^ 



^72 3V3 ^cTf - ^v/2 3^/^ hi 



MO) 



33 



-10 



g^'-' - 18^2 ^0 



fe(i) 



-^v/2 



.(1) 

-10 



,(0) 



+ a« + 24v^ af ) + ^^v^ h^^ - 12v^ /.« , 



(2.67) 



16 



3456 ^ ° ^ 9 ^^ 1296 ° 16 ^° 384 ^° 

16 144 Z 

(2.68) 

Given data in Table 1 we can verify (2.63)-(2.64) to agree with a relative accuracy 
of ~ 10"^°, and (2.65)-(2.67) with a relative accuracy of ~ 10^'^. 

We are now in position to evaluate the free energy near the critical point. First of 
all, notice that i^oT^ or ^^ will receive 0{\^) corrections: 

, _o , ^' ^512v^ ,(0) 320v^ ,(0) 184^2 ,/6 .(o) 8^2 ./^ / mo)^ 

112V^ /6 .HI) , ^^ Ml) , 1024V6 ,(1) 128v^ ,(i) 
"1 ^^ "J '^0 "1 ^^ "2 "1 ^^^ y2 5 yo 



i^3v«cr'er') H- o (m; 



TT \ yU^ \9 27 3 6 



(2.69) 



(2.70) 



(2.71) 

with the mass deformation lowering^ the critical temperature for a fixed chemical 
potential. Once again, introducing 

t=--l, /i = l-^, (2.72) 

Tc He 



^We used results of Table 1. 
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we find 



(2.74) 



For the first two corrections to (2.49) we have 

1 1 



3? = - 2 Af - 5 + ^V6 cS> , (2.75) 

,; . _ 8V2 /4" + ?V2 M»' - i| 3V« (cj'"')' - I 3=/« 4"»cJ'" + H!i 3'/« aJW 

+ f 4""-|v2 + iv3ci;', 

(2.76) 



A(0) 



.» = - 4 M" + iV6eS' - ^ - f V6 gS'»' + ^V6 gj'»' + | 3'/« V2 <' 

+ ^V2ar>-?3V«V2(cr)', ''''" 

si = - 16^2 M" + f 72 /4»' - ^ 3V« (4'"')' - I 3=/«' cr'4'" + ^ 3'/» ftj'" 
736 1/8 m 512^ <.(i) , 2048^ ,.,„ 4096^ j,r, 1280 ,.,„, 



81 ^ 27 "^ 243 "^ 81 "" 27 

(2.78) 

It is straightforward to clieck, tliat given constraints from tlie basic tliermodynamic 
relations (2.63)-(2.67) 

wliicli guarantees tliat tlie mass-deformed tlieory lias the same static critical exponents 
(a,/?, 7, (5) as in (2.53). 



1^ 



3 Fluctuations in charged plasma and its holographic dual 

3.1 Hydrodynamic modes in charged plasma 

In this section we consider the propagation of hydrodynamic modes in four- dimensional 
charged relativistic plasma in the absence of background electro-magnetic fields. 

The first order hydrodynamic equations of motion in Minkowski space-time are 
simply the conservation laws for the stress-energy tensor and the f/(l) current: 

d^T^"' = 0, 9^J^ = 0. (3.1) 

One can do the standard decomposition of the stress tensor, 

T^"" = eu^v' + PA^'^ + W , (3.2) 

where 

^f^i^ = ^M'^ + u^'u'' , WX = , (3.3) 

and e and P are the energy density and the pressure respectively. The dissipative term 
YiiJ-i^ is given by 

n^'^ = -ria>'' - CA^'^(9«m") , (3.4) 

where 

^t^u = A^"A'^^(9,M/3 + dpda) - ^A'^^A'^^idX) , (3.5) 

and 1] and ( are the shear and the bulk viscosities. Note that 11'^'^ is, by definition, zero 
at local equilibrium. The current J^ is given by 

J^ = pu^' + pt" , (3.6) 

where v^ is the dissipative part satisfying ■u'^z/^ = 0: 

y>^ = aQA^''{~^,^^+^^,T) . (3.7) 

In this expression, T is the temperature, /i is the chemical potential, and aq is the 
conductivity coefficient. We would like to study fluctuations around the equilibrium 
state in which 

«'' = (1,0,0,0), T = const. , /i = const. . (3.8) 

As an independent set of variables we will choose the three spatial components of 
the velocity 5u^ = 6u^, 5v? = Su^, 6u^ = 6u^, as well as 6T and 6fj,. As usual, all 
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perturbations are of the plane-wave form exp{—iut + iqz). We find that the relevant 
fluctuations of T^'^ are 

ST'' = (e + P)5u' , 
Similarly, we obtain the following fluctuations of the current 



(3.9) 



,T" = 5P-{-,^C) dM' = yj-^^i, + [^-^j^ST - ^-, + C ) d,iu. 



5.r = p5u' + aq {-djp + ^d,5T 



(3.10) 



Substituting these expressions into equations of motion (3.1) and performing a Fourier 
transformation we get the following system of equations 

0=u(t + P)5n'-qU^\ 6^+(^^\ dT]+iq'(\r, + OSu', (3.11) 



=^ ■ • A 



^^ + {%) ^A - IP^""' + 'I'^'Q {^P - J^^^ 



0=fa;(e + P) + ig2^j(5M% 
0=(w(e + P) + ig2^ )5M^ 



(3.12) 



where we assembled equations in the decoupled sets. The three equations (3.11) de- 
scribe the propagation of sound waves in charged plasma, while the pair (3.12) describes 
the propagation of the shear modes of different polarizations. 
In the shear channel, (3.12), the dispersion relation takes form 

uJ = tqV = -tq ——- = -1—-- . 3.13) 

t + F I s 1 s — pp 
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Notice that even though in a holographic gauge theory plasma we study here the ratio 
of the shear viscosity to the entropy density is universal [19] 

s An 
a combination involving the diffusive constant TV has a nontrivial dependence on - 
and ^. 

In the sound channel, (3.11), the dispersion relation takes form 

u = ±c,q - iTq^ + 0{q^) , (3.15) 

where the speed of sound Cs and the attenuation F are given by 



c 



diP,p) , 9(e,P)\/, , „, die,p) 



d(T,p) ' diT,fi)jy ' d{T,p) 



e + P)-)-l^ + p^^)-H (e + ^)T^7FT > (3-16) 



-1 



-1 






dlnfijrp \d\nT J I \\d\nfj,Jj, \d\nT ^ 

(3.17) 

It is instructive to analyze (3.15) for conformal theories and for theories with softly 

broken scale invariance, i.e., when 

MM , , 

— «1, — «1, (3.18) 

where M is a conformal symmetry breaking scale. 
■ In a conformal hydrodynamics 

C = 0, e = 3P, e^T^^(l), p ^ T=^ 7^ (^) , (3.19) 

where S and TZ are functions of the dimensionless ratio ^, thus 

d{P,p) _1 d{e,p) d{e,P) _ 1 die,e) _ 

d{T, p) 3 (9(T, p) ' a(T, p) 3 (9(T, p) 

(5^)/(atV)„--'^-'--'---'(-?^)--^- 



•^O , / -96 ^^ _y4^^,^4y4^^y5 / ^1 ^/^4^^ 



ain/xy^ V-^inTy^ t v t^ 

aP\ ^ .9(P,e) ^ 1 %,e) ^ Q 
dp J , d{p,e) 3 9(p,e) 
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(3.20) 



As a result, we find for tlie speed of sound and tlie sound wave attenuation 

CFT "^ ' CFT 3(e + P) 3T s Ts — jip 

Once again, notice that TV lias a nontrivial dependence on ^, even though the shear 
viscosity ratio is universal (3.14)^. 
■ When the scale invariance is softly broken (3.18), we find 



2 1 .^ (M^ 



' 'l + — +-x-^xC — - xC — - 



(3.22) 



3(e + P) V 4:7] J T T \T^ J \T^ 

where, given (3.20), we used 

In (3.22) and (3.23) we suppressed dependence on ^ in 0{M'^) terms. What is impor- 
tant for the later discussion is that the contribution of the term oc cjq to the attenuation 
risoforder C(M4). 

3.2 Sound waves of mass-deformed RN plasma 

We can use (3.16) and (3.17) to compute the speed of sound and the bulk viscosity of 
the mass-deformed RN plasma (2.1) to order 0{\^). 
We find it convenient to introduce 

"-if- '-if' P-^*' 

where {cj, g' = |^} are the frequency and the momentum of the sound mode. Further- 
more, we parametrize the sound wave dispersion relation (3.15) as 



ra = -^ /3i - ^^ /32 + O(q') 






(3.25) 



n=0 ^ ^ n=0 



^In conformal holographic hydrodynamics this was first estabhshcd in [20] 

22 



In this parametrization /3i = /32 = 1 for a conformal plasma with vanishing chemical 
potential [21]. Using (2.42), (2.43), (2.59), (3.14) and introducing 

/3 



ao 



yrr 



K 



ai = -/3V2 + \'^(3 ai,2(/t), 



'.^^^^h.J.^ n. = R r:=\ (3-26) 



/ii = 7-7== + A^ /ii,2(/t), 90 = P, Co = A, 
3vl + /t 



cio = ciq{k) , 
for the UV parameters (2.30), and 

4 = ^f, r-— + A^/3 4M , (?o' = /3 + A'/3 ^^l^:) 
92 = A'/3 ^2'2('^) , h'o = (1 + '^)'/' + A' <2('^) , 

for the IR parameters (2.37), we find from (3.15)-(3.17) 

K + 2 

/5i,o = ±1 , /32,o = 7^ — -77 , 

ZK + 2 

fiV2 1 



(3.27) 



(3.28) 



while /32,i directly determines the bulk viscosity to the shear viscosity ratio, since 
contribution to the attenuation proportional to the conductivity aq vanishes to order 
(9(A^), see (3.22), and the conductivity itself is finite at criticality [7], 

where 

,, , 8(1 + k) ^ 4(l + /t)l/22l/2 4 xl/2 /, . X 

Al'^) = ^^^^ /32,i - ^ ^ /^i,2(«:) + 3(«: + 2)(1 + /.)i/2 ^^_^(^) 

20 , ,,^, , 4(«: + 2)2 , ^^^^ , 2(5/. + 4) 



6 6(1 + K 



(1 + K]^i'^K 

(3.30) 



In the next section we study sound channel quasinormal modes of mass-deformed 
RN black holes (2.5) and compute (5i and /32 in (3.25) to order 0{\^) inclusive, at 



K = 2. Since cio 



= c\q in Table 1, (3.28) would provide a highly nontrivial test on 

K = 2 
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a consistency of our thermodynamic and hydrodynamic analysis. Moreover, given that 



hi: 



92,2 



K = 2 



K=2 



'n 5 



9^ 



'-2,2 



MO) 



K=2 



'*0,2 



h 



h{0) 
) 



^2,0 



"0,1 



K = 2 



K=2 



K = 2 



.;^°^ 



>(0) 
-0 ; 



(3.3i: 



as presented in Table 1 are finite, any possible divergence in the bulk viscosity at 
criticality, i.e., aX k, = Kc (2.69), would arise from the divergence of /32,i = P2,i{i^), 
as K — 7> 2. The bulk viscosity to the shear viscosity ratio to leading order in ^r at 
criticality is then computed from (3.29): 



lim - = A^ lim AU) + C(A^) 



^^-hmA(.) + 0U^ 



(3.32) 



3.3 Fluctuations of the deformed RN black hole 



To determine dispersion relation of the sound channel quasinormal mode we have to 
analyze fiuctuation in the background geometry {g^u,Afi,(p,x} (2-5): 



M ' 



A^-^A^ + 6 A 
(f)^ (f) + 6(f), 



(3.33) 



We choose the gauge 



h 



tr 



h 



XiT 



hr 



0, 



5A, = Q. 



(3.34) 



Additionally, we assume that all the fiuctuations depend only on (t, x^, r), i.e., we have 
an 0(2) rotational symmetry in the a;ia;2 plane. At a linearized level, the fiuctuations 
{hfj^y, 5A^, 6(f), 6x} of different helicities with respect to this symmetry will decouple 
from each other. The sound channel quasinormal mode corresponds to helicity-zero 
fiuctuations [22]: 



{htt ; haa = "-xx + hyy , htx3 , hx^xa 5 oAf , oAx^ , b(p , bxs 



(3.35) 
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We introduce 

htt =(r^'^''^^ c\ Htt , 

h —p-i'^t+iqx'i 2 TT 
'Hx3 — ^ ^2 tz ■} 

U —p-ii^t+iqx;i 2 tt 
"'aa ^ "^2 ^^aa i 

h _ -iut+iqxz 2 TT 

"'xyx'i — c 1^2 ^-^X'iX-i , 

(3. do) 

KA —p-ii^t+iqxz A 

where {Htt,Htx.^,Haa,Hx.^x3,At,Ax3,p,c} are functions of a radial coordinate only. 
From the effective action (2.1) it is straightforward to derive 8 second order differ- 
ential equations of motion for the fluctuations, and 4 flrst order differential constraints 
associated with fixing the gauge invariance as in (3.34). Altogether we expect 8 — 4 = 4 
independent gauge-invariant combinations of fluctuations. Indeed, analyzing the trans- 
formations of (3.36) under the residual gauge and diffeomorphism transformations it 
is straightforward to construct these combinations: 

/~f in c /^ \ n c 

'Z'H =4— Htz + 2 Hzz — Haa 1 7^—, — + 2^^ Htt , 

CO \ UJ"^ C2C2 / CO C2 



q [In 4 



Za —-^t + — •A.xz — T, 4T7 Haa 5 



1/ 



(3.37) 



^^-P^^" 



^ Ha 



[Inct]' 

With somewhat tedious analysis we can verify that equations of motion for the gauge- 
invariant fluctuations {Zh, Z^, Zp, Z^} decouple. 

The spectrum of quasinormal modes is determined [22] by imposing on {Z^^ Z^-, Zp, 
Zc} an incoming wave boundary condition at the horizon, and requiring vanishing of 
the non-normalizable modes for {Zh, Z^, Zp, Zc} near the boundary. In the hydrody- 
namic limit tt) — )■ 0, q — ?■ with — kept flxed, this leads to the following perturbative 
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expansions 



^H 



[l-x) 



ZHfl + n zh,i + 0{cC 



Zj, =(1 - x)-*" ^^,0 + ^q ZA,i + 0(q' 



^. =(1 - x) 



:i-x) 



Zpfl + i(\ Zp^i + C(q^ 



Zcfi + «q 2c,i + C>(q^ 



(3.38) 



with the following boundary conditions on {zH^i,Za^i.,Zy^^i\: 



lim ZHfl = 1 , lim zh,i = , lim zj^^i = lim Zp^j = lim Zc,i = finite , 

r — ^1 rr — ^1 rr— ^1 .t — ^1 .t— ^1 



^H,^ = 0(x) 



X^l- 



X— s-l- 



X— S-l- 



^A. = O (xV2) , ^^^, = O {x'/') , z,, = O (x^/^) , 



;e- 



invariant 



(3.39) 

as X — )■ 0+. 

To leading order in the hydrodynamic approximation, wave functions of the gaugi 
fiuctuations {zha, zj^fi, Zpfl, Zcfl} satisfy the following equations 

~^Hfl + ^101 ZuQ + Ci02 Zji^^Q + Ci03 ^p,o + ^104 ^c,0 + ^105 ZHfl + Ci06 ^^,0 

+ Ci07 ^p,0 + Ci08 ^c,0 5 
~^yl,0 + ^201 Z^Q + C202 ^^,0 + ^203 ^p,o + ^204 ^c,0 + ^205 ZHfl + C206 Z_Afl 

+ C207 ^p,0 + C2O8 ^c,0 5 
~^p,0 + ^301 Z^Q + C302 Zj^Q + C303 ^p + ^304 ^c,0 + ^305 ^H,0 + C306 ^yl,0 

+ C307 Zpfl + C308 ^c,0 5 
~^c,0 + ^401 ^H,o + ^402 Zji^Q + C403 2p + C404 ^c,0 + ^405 ZHfl + C406 ^^,0 

+ C407 Zpfl + C408 Zcfl , 

where connection coefficients Cjoj are complicated nonlinear functionals of the back- 
ground fields {(7, A, 0, x} with explicit dependence on x and /3i2 = /^^ , and implicit 
dependence on k and A: 



(3.40) 



^iOj ^iOj 



{fi',^,0,x}; a:; /3i2 



(3.41) 



To leading order in the hydrodynamic approximation, wave functions of the gauge- 
invariant fiuctuations {zh,i, Zji,^i, Zp^i, Zc,i} satisfy equations identical to (3.40), apart 
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from the source terms {Jh, Ja-, Jpi Jc}'- 



Cioi Zjj^ + C102 Zji^^i + C103 2p,i + Ci04 ^cl + ^105 Zh,1 + C106 2yl,l 



— ^// 1 + Cioi Zii^^ 

+ Ci07 ^p,l + C108 ^c,l + iTh 



=Zj^^x + C20I ^j:/,! + C202 ^^,1 + C203 ^p,l + ^204 ^c,l + ^205 ^i/,1 + ^206 ^^4,1 

+ C207 ^p,l + C2O8 ^c,l + J A ■, 
~^p,l + ^301 ^H,l + ^302 ^yl,i + C303 ^p 1 + C304 ^c,l + ^305 ^H,l + C306 ^yl,l 

+ C307 ^p,l + C308 ^c,l + i^ 5 
=^c,l + ^401 Zjji + C402 2_4^i + C403 2p 1 + C404 ^cl + ^405 ^i/,1 + C406 ^yl,l 

+ C407 2p,l + C408 ^c,l + >Jc , 











(3.42) 



with 



>Jh =Ciii ^ji^o + C112 2_4 + C113 2;^ Q + C114 Z^ Q + C115 Zfjfi + C116 Z_4^o 

+ Cm Zpfi + C118 ^c,o 5 

1Z4 =^211 2//^o + ^212 Zj^Q + C213 2p + C214 ^cO + ^215 Zff^Q + C2I6 2^4,0 

+ C217 Zpfi + C218 ^c,0 5 
Jp =^311 2j^o + ^312 ^^,0 + ^313 ^p,0 + ^314 ^c,0 + ^315 ^Hfi + C316 ^^,0 
+ C317 Zr,,o + C318 ^c.O 1 



(3.43) 



+ C317 2;p^o + C318 Zcfl 
Jc =Ca,\\ Zjjq + C412 2_4 + C413 z^^Q + C414 2;^ Q + C415 2;//^o + C416 2^,0 
+ C417 2;p^o + C418 Zcfi . 

The new connection coefficients Cnj are comphcated nonhnear functionals of the back- 
ground fields {g,A,(f),x} with exphcit dependence on x and P12, f32, and imphcit de- 
pendence on K, and A: 



c 



ilj 



Wlj 



{g,A<P,x}; x; {/3i2,/32} 



(3.44) 



Exphcit expressions for Cikj are available from the author upon the request. 

It is straightforward to determine the asymptotic expansions (satisfying (3.39)) 
and thus set up the boundary value problem what would determine {/3i2,/32}, along 
with the coefficients of the normalizable modes for {zH,i, ZA,i, Zp^i, Zc,i}- We found that, 
unfortunately, directly solving the resulting boundary value problem (as it was done 
for example in [23]) is not possible with reasonable computational resources. In the 
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rest of this section we explain the origin of the problem and outline the solution^. 

Recall that (3.40) and (3.42) were obtained in the hydrodynamic limit from the 
full quasinormal mode equations for {Zh, Zj{, Zp, Zc}. The latter equations have a 
structure identical to that of (3.40). For example, we have 

=Zfj + Ch,i Zfj + Ch,2 Z^ + Ch,3 Z + Ch,4 Z^ + Ch,5 Zh + Ch,6 Zj{ 

(3-45) 

+ Ch,7 Zp + Ch,8 Zc , 



with 



Cha — C 



HA — ^HA 



{^,A0,X}; x; {tTJ,q} 



(3.46) 



Explicit expressions for CH,i show that (some of) these coefficients have a simple pole 
at 

= tt)2-(l-a;)2q2. (3.47) 

The residues of these poles are always proportional to the derivative of the background 
gauge potential^ A', or {A'^, (see (2.5)), and does not vanish when either (or both) 
= 1 and X = 0. In other words, these poles always occur in studies of the sound 
channel quasinormal modes of charged black holes in asymptotic AdS^ geometry. Since 
the speed of sound squared c^ in holographic plasma with a UV conformal fixed point 
is bounded by [26, 27]^ 

4<l, (3.48) 

from (3.47) we always expect to encounter a singularity in the connection coefficients of 
the differential equations describing the propagation of the sound channel quasinormal 
modes of the charged black holes inside the range of the integration, i.e., for x G 
(0, 1). Of course, the singularity in the coefficients of the differential equation does 
not imply that the solutions are singular — in our case we find that they are not — 
however, the presence of such singularities poses technical difficulties for the boundary 
value problem one has to solve. Actually, the issue is even more complicated: the 
successive hydrodynamic approximations produce equations ((3.40) and (3.42) in our 
case) which connection coefficients involve successive derivatives with respect to Cg. 
Thus, while some of the coefficients Cjoj in (3.40) have a simple pole inside the range 



^A simple explicit example explaining the use of a new computational technique is presented in 
Appendix B. 

*^This explains why such poles are absent in the analysis of [24] or [25]. 
^See [28] for an exception. 
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of integration, some of the coefficients Cnj in (3.43) have a double pole! Finally, a 
perturbative expansion in the deformation parameter A, which is needed to extract 
/32,i and ultimately the bulk viscosity of the charged plasma, see (3.32), produces yet 
higher order poles — here we need to deal with the third order poles in the connection 
coefficients inside the range of the integration in the boundary value problem. We find 
that numerical techniques for solving the boundary value problem developed in [29] 
become unreliable once the connection coefficients of the differential equations have a 
second or higher order poles. A new approach is needed. 

Basically, we need to reformulate the boundary value problem in such a way that 
connection coefficients of the corresponding differential equations have simple poles 
inside the range of integration, at worst. First, even though we are after the hydro- 
dynamics to order (9(A^) only, we treat (3.40) and (3.42) exactly in A. The boundary 
value problem would then determine 

/9i = /3i(A), /32 = /32(A). (3.49) 

From the data sets for small A we can extract /3i,05/3i,i, /32,o? /32,i (see (3.25)). Such a 
step insures that Qqj have simple poles at most, while Cnj have double poles at most. 
In particular, the boundary value problem for (3.40) is amenable to the treatment of 
[29]. Second, the hydrodynamic origin of the second order poles in Cnj implies that 
the residues of the second order poles in Cnj and ^^ must be proportional to each 
other. Indeed, for each value of indexes {i,j} we find 



%~2 d 



Ciij - -2/32 /- 'Q0~^iOj + Ciij , Ciij - Ciij 



{g,A(t>,x}; x; /3i2 



(3.50) 



with Ciij having at most simple poles for x G (0, 1). Given (3.50), the boundary value 
problem (3.42) can be reformulated as follows. 
■ We represent 



zh,i = - 2/32 

ZA,i = - 2/32 

(3.51) 

Zp,l = — 2/32- 



ZcA = - 2(3: 



V/3l2 


d 

5/3l2 


zh,c 


+ Zh,1 


V/3l2 


d 

d/3i2 


ZA,(] 


+ ZA,l 


V/3l2 

V3 


d 

d/3i2 


Zpfl 


+ Zp,l ) 


V/3l2 


d 

d/3i2 


Zcfi 


+ ~Zc,l ■ 


2c 


) 







■ The wave functions {zh,0: ^a,o, ^p,o, ^cfl} satisfy exactly the same equations as {zh,o, 
Zj\,^o,Zpfi,Zcfi}, i.e., (3.40), with f3i2 being treated as an extra free parameter and the 
only change in the boundary conditions being 

hm ZHfl = finite = 2o(/3i2) • (3.52) 

x— i>0+ 

Clearly, imposing the Dirichlet boundary condition Zq{(5i2) = would determine the 
value of /3i2 which would identify all z...^q with z...^q: 

Z-,0 =z...fi, (3.53) 

Zo(l3i2)=0 

where ■ ■ ■ stand for any of {H,A,p,c}. For each value of A, we can evaluate Zq for 
some set of /3i2 and compute 



d 



(3.54) 

2o(/3i2)=0 



The boundary value problem(s) implementing this procedure involve solving ODE's 
with simple poles at most in the connection coefficients inside the range of integration. 

■ Given (3.51) and the definition of z...^q, it is straightforward to see that z...^i would 
satisfy equations identical to (3.42), except that now the source terms (3.43) would be 
constructed from the connection coefficients Cnj. For example, 

>Jh -^ Jh = Cm ZffQ + Cii2 Zj^Q + Cii3 z Q + Cii4 z^q + Cus zh^o + Cue Zji^^Q 

(3.55) 

+ Cii7 Zp^o + Cii8 Zc,0 . 

By construction, this final boundary value problem would involve ODE's with simple 
poles at most in the connection coefficients inside the range of integration. However, 
since it does not depend on /32 (see (3.50)), to find a solution we must modify^° the 
Dirichlet condition at the boundary for zh,i'- 

lim Zhi = finite = Zi . (3.56) 

x— )-0+ ' 

■ Ultimately, zh,i must satisfy the Dirichlet condition at the boundary (3.39). Using 
(3.51), (3.52), (3.54) and (3.56) this leads to 

-2/32^^^ + ^1 = 0, (3.57) 



''This is the only modification in the boundary conditions. 
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which determines /32 as 

/32 = ^ § . (3.58) 

Thus, we succeeded in reformulating the boundary value problem for computing /32 in 
such a way that it does not involve numerical integration of ODE's with connection 
coefficients having singularities stronger than simple poles for x G (0, 1). 

We emphasize once again that the procedure described above is necessary only 
in computing sound channel quasinormal modes for charged black holes in asymp- 
totic AdSr, geometry. On the other hand, it is generic, and can be applied to AdS- 
Schwarzschild black holes as well. In Appendix B we demonstrate the new method 
for the computation of the dispersion relation of the sound waves in strongly coupled 
A/" = 4 SYM plasma. 

3.4 Mass-deformed RN plasma transport at criticality 

In this section we present results for the speed of sound and the bulk viscosity of 
the mass-deformed RN plasma in the vicinity of the second order phase transition. We 
perform numerical analysis for a fixed /t = 2 as a function of (a small) mass- deformation 
parameter A, see (2.59). In general (see (2.69)) 

Kc = Kci\) = 2 + 0{\'^). (3.59) 

However, setting /t = 2 is still sufficient to extract the leading correction to the speed 
of sound 

^hm^ (cl) + O(A^) = ^ mn = 2, A))^ , (3.60) 

and the leading contribution to the bulk viscosity near the vicinity of the phase tran- 
sition, see (3.32): 



lim A"' X - +C(A') = lim A(fi; = 2,A). (3.61) 

3.4-1 The speed of sound 

First, we verify that in the conformal limit the speed of sound is independent of the 
chemical potential, see (3.21) . We find 



mK,\ = o)f-i 



10~'^---10~'S KG [0.4,2]. (3.62) 
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Figure 1: (Colour online) Deviation of the speed of sound (1 — f3f), see (3.25), in mass- 
deformed RN plasma from its conformal value as a function of the mass- deformation 
parameter A at k = 2, see (2.59). The blue dots represents results obtained from the 
holographic hydrodynamic equations (3.40), and the solid red line represents thermo- 
dynamic prediction, see (3.28). 



Figure 1 presents the results for the deviation of the speed of sound in RN plasma 
from the conformal value as a function of the mass-deformation parameter A: the blue 
dots are obtained from directly solving the holographic hydrodynamic equations to 
leading order (3.40), while the solid red line^^ 



(1 - /3i ) 



-2/3i,i A^ 



0.061132('8) 



(3.63) 



red 



is the thermodynamic prediction for this deviation, valid in the limit A — ?■ 0. The 
results are in excellent agreement: for instance, for A = 0.01 we find 



:i - /3i 



red 



y^ Pi) blue 



-1 



5 X 10 



-6 



(3.64) 



3.4-2 The bulk viscosity 

Scale invariance of the model at A = predicts (3.28) 

2 

/32(A = Ojprediction = 1^2 fi\^^ ~ '^' ~ '^ ' 



(3.65) 



^We used (3.28) and the result of Table 1: cio(k = 2) = 



-0\- JO) 
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Figure 2: (Colour online) Blue dots represent Zq{I3i2) (see (3.52)) for the mass- 
deformation parameter A = 0.01. /3^2 — (/5i('^ = 0.01)) . The solid red line represents 
the best linear fit to Zq. 

Explicit computation of the attenuation coefficient P2 following the method explained 

in section 3.3 yields 

/32(A = 0) 



. ,- 2xl0-^ (3.66) 

y2\^ ^Jprediction 

Before we present a general plot for the bulk viscosity in RN plasma, we discuss in 
some details results for A = 0.01. Blue dots in Figure 2 represent the results for Zo(/3i2) 
(see (3.52)) for values of /3i2 in the vicinity of /3*2 = (/3i(A = 0.01)) as evaluated from 
solving the hydrodynamic equations at leading order (3.40). We present the data as a 
function of (/3i2 — /3i2) since, following the discussion around (3.52), we expect 



2o(/3i2 = (3 



12 J 



0. 



(3.67) 



The solid red line represents the best linear fit to the hydrodynamic Zq data. Explicitly, 
we find 



Z, 



8. X 10"'^ - 1.400011(9) (/3i2 - 13*2] 



red 



From (3.68) we determine 



7' 



-1.400011903986774, 



A=0.01 



Next, we evaluate Zi as explained in section 3.3: 



Zi 



-1.0777132662397266, 



(3.68) 



(3.69) 



(3.70) 



A=0.01 
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Figure 3: (Colour online) Blue dots represent /32(A) (see (3.58)) at k = 2. The solid 
red line represents the best quadratic fit to the first 10 blue dots. 



which following (3.58) computes ^2'- 



^2 



0.6666585596689777. 



(3.71) 



A=0.01 



Figure 3 presents results (blue points) for /32(A), were we simply iterated the proce- 
dure described above for /32(A = 0.01). The solid red line represents the best quadratic 
fit (with the linear term absent) to the first 10 blue points: 



/32 



0.666666(3) - 0.081002(5) A^ . 



(3.72) 



red 



Note that /32(A = 0) in this fit is remarkably close to exact CFT value (3.65). We can 
now compute /32,i in (3.25) as 



/3: 



1 rf2 



2,1 



2d\'^ 



(3, 



/32,i = -0.081002(5) . 



(3.73) 



red 



Finally, using results in Table 1, we find from (3.30) 



thus, from (3.29) 



A = 0.224133(3) + 2/32,i 



0.062128(2) A2 + C(A^ 



(3.74) 



(3.75) 
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Curiously, since 

\-cl^ ^(-2/3i,i) = ^(1 - /3i^2) = 0.020377(6) A^ + 0{\') , (3.76) 

the bulk viscosity bound for holographic gauge theory plasma proposed in [13] 



^>2^3-e,^j, (3.77) 

is satisfied in the vicinity of the phase transition of mass-deformed RN plasma. 
To summarize, we find: 

^ = 3.0488(5) (^ - ^') + ^ ( (^ - ^') j ' (3.78) 

for the bulk viscosity of the mass-deformed RN plasma in the vicinity of the second 
order phase transition. 



4 Dynamical critical phenomena in RN plasma 

In this section we study the response of conformal and mass-deformed RN plasma to 
inhomogeneous and time-dependent variation of the chemical potential 

At a linearized level the variation of the chemical potential would produce a correspond- 
ing variation in the charge density, 6p{t,x) ( p^^ for the Fourier components). As in 
conventional theory of dynamical critical phenomena [8], we introduce the dynamical 
susceptibility x^ k, 



X.i. 



The fiuctuation-dissipation theorem states that 






(4.2) 



Giu,k) = — lmx^,:, (4.3) 

UJ 

where G{uj, k) is a Fourier transform of the charge density variation two-point correla- 
tion function 

G{t,x) = {6pit,x)6p{0,0))^=o. (4.4) 
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Furthermore, the equal-time correlation function 

G{k) = G{u = 0,k), (4.5) 

is related to the static susceptibility 

Xk = X^=o,k, (4.6) 

by the equipartition theorem 

Gik)=Tx^,. (4.7) 

In the vicinity of (but not at) the critical point t = ^ — 1— 7'0(t7^0) the equal-time 
correlation function G{x) decays exponentially 

G{x) oc e-l^^l/« , (4.8) 

where ^ is the correlation length, implying that G{k), and through the equipartition 
relation (4.7) the static susceptibility x^, have a pole at 

P oc -r' • (4.9) 

Right at the critical point, t = 0, the equal-time correlation function has a power-law 
decay 

G'(f) oc |f 1"^+'-^' ^ Gik) oc \k\-^^\ (4.10) 

where p = 3 is the number of spatial dimensions and rj is the anomalous scaling 
exponent. 

The theory of dynamical critical phenomena [8] predicts that in the vicinity of the 
continuous phase transition and for \k\ ~ C,^^ the full dynamical susceptibility x^i k ^^^^ 
develop a pole at 

a;oc-ir% (4.11) 

with z being the dynamical critical exponent of the system. The frequency in (4.11) 
defines a relaxation time r~^ 

T~^=iuoci-\ (4.12) 

characterizing the equilibration time scale of the dynamical system. 

In the rest of this section we analyze dynamical susceptibility of the strongly coupled 
conformal RN plasma. Following up the poles in static susceptibility (see (4.9)) in the 
vicinity of the phase transition we determine the (static) scaling exponent v of the 
correlation length: ^ oc t^'^ . The scaling of the pole in the dynamical susceptibility 
(see (4.11)) determines the dynamical critical exponent of the RN plasma. Finally, we 
comment on dynamical critical phenomena in mass-deformed RN plasma. 
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4.1 Computation of x^j: in holographic dual 

Without the loss of generahty we can assume that 

k' = q6i. (4.13) 

In dual gravitational description the variation of the chemical potential fi^^^k trans- 
lates into the variation of the non-normalizable mode of the fluctuation of the bulk 
vector fleld At (3.36), correspondingly the variation of the non-normalizable mode in 
the gauge-invariant fluctuation Z^ (3.37). Up to an overall factor^^ (see (2.32)), the 
variation in the charge density p^^^k is the normalizable component of At (or Z^ for 
the gauge-invariant fluctuation). The fluctuations of the dual gravitational background 
were analyzed extensively in section 3. As we study here the dynamical critical phe- 
nomena of the conformal RN plasma, we can consistently set Z^ = 0. We are left with 
the linear coupled system of fluctuations {Zh, Z_a, Zp} (3.37). There are two important 
differences in analysis of this system of fluctuations compare to the one in section 3, 
relevant for the computation of the sound wave dispersion relation: 

■ In both cases {Zh, ^^, Zp} must satisfy an incoming boundary conditions at the 
horizon. In case of the sound waves, we had to impose the vanishing of all the non- 
normalizable modes for {Z^, ^^, Zp} at the boundary (see (3.39)). As a result, solving 
the boundary value problem determined the dispersion relation for the sound waves 
(3.25) 

m = m(q). (4.14) 

In computing the dynamical susceptibility, both vo and q are independent; nonetheless, 
the boundary value problem on the gauge-invariant fluctuations has a solution because 
the non-normalizable component of Z_a at the boundary is nonzero now. 

■ Transport coeflicients are encoded in the perturbative in q expansion of the sound 
wave dispersion relation (3.15). Thus, it was suflicient to implement the perturbative 
(hydrodynamic) expansion for the fluctuations (3.38) and the dispersion relation (3.25). 
In study of the dynamical critical phenomena of strongly coupled RN plasma we are 
interested in the poles (for complex tt) , q) of the dynamical susceptibility Xm,q- Thus, 
we can not do the computations perturbatively in VD,q. Indeed, since the system 
of equations of motion for the fluctuations {Z^, Zj^, Zp} is linear, without the loss 
of generality we can set the non-normalizable component of Z^ near the boundary. 



"'^^This factor is finite in tlie vicinity of the transition. 
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namely /i,^_q to one: 

/X„,q = l. (4.15) 

From (4.2) we have 

Xw,q=Proq, (4.16) 

thus any poles in the dynamical susceptibility must come from the poles in the normal- 
izable component of Zj, near the boundary, and as such they must be non-perturbative 
in m , q. 

We can summarize now the boundary value problem whose solution would deter- 
mine the dynamical susceptibility. Introducing^'^ 

Zn =(1 - x)-''^ VQ-^ Zh{x, It), q) , 

Z^=il-xr"'q~'zAx,w,q), (4.17) 

Zp={l-x)-''° q~2 Zp{x,K>,q), 

the equations of motion for {z^, z^, Zp} are solved with the following boundary condi- 
tions: 

lim Zh = lim 2;^ = lim Zp = finite , 

zh = 0{x), Zj^ = l + Z{rv,q) x^^^ + 0{x), Zp = 0[x^l^), as a; ^ 0+ . 

(4.18) 

The normalizable component of Zj^ near the boundary Z is proportional to the dynam- 
ical susceptibility: 

X„,q oc Z(tt), q) . (4.19) 

4.2 Static and dynamical susceptibilities of RN plasma 

As discussed in section 2.5, a second order phase transition in strongly coupled RN 
plasma happens for 

K = 2 ^^ ^ = ll, (4.20) 

/i vr 

We do not provide the technical details of the analysis of the boundary value problem 

(4. 17)- (4. 19) and present only the results. We point out that the coefficients of the 

differential equations for {^z^,^ Zj^^ Zp\ are non-singular as k — )■ 2; thus the appearance 

of poles in Z in this limit is not obvious. 



^■^Thc ro— and q— dependent rescaling are for convenience in further analysis. 
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Figure 4: (Colour online) The scaling (blue dots) of the inverse of the static suscepti- 
bility Xro=o.q=o in the vicinity of the critical point. The solid red line is a quadratic fit 
to the data. 

Figure 4 shows the inverse of the static susceptibility at q = (blue dots) in the 
vicinity of the critical point. The solid red line represents the best quadratic fit to the 
data: 

Z^.j^ = 4.74744 ■ 10^^ + 0.125116 (ft -2) -0.0156647 (k- 2)2 + C((/t- 2)^). (4.21) 

The red line (4.21) intersects the k axis at 

Kc = 1.999999(6) , (4.22) 

in excellent agreement with the expected value Kc = 2. Thus we reproduce the ther- 
modynamic result for the static susceptibility 



Xw=0,q=0 = XT OC Z OC 



1 



OC +r 



-1/2 



rL Kj(^ ^^^^ rxiQ 



(4.23) 



where we used relation (2.50) between k, and the reduced temperature t. 

Figure 5 presents the poles (blue dots) of the static susceptibility at q = q* in the 
vicinity of the critical point: 

X^io,q=q, = . (4.24) 

The solid red line represents the best quadratic fit to the data: 

g2^.^ = -2.32509 • 10^^ + 0.11873 (k- 2) -0.0347648 (k- 2)2 + C((/t- 2)2). (4.25) 

39 



0.04 



0.02 



-0.02 



-0.04 




Figure 5: (Colour online) Poles of the static susceptibility in the vicinity of the critical 



point: X 



-1 

rD=0,q=q* 



0. The solid red line is a quadratic fit to the data. 



Notice that in the stable phase, i.e., for k < 2, in the vicinity of the phase transition 
the poles in the static susceptibility are for purely imaginary momenta, which im- 
plies the exponential decay of the charge density two-point correlation function (4.8). 
Furthermore, from (4.9) we identify the correlation length as 



(27rreO' oc q;2 oc 



1 



K Kc 



OC -f f 



-1/2 



Q < Kc — K <^ Kc 



(4.26) 



where we used the results of the fit (4.25) and the relation between k and the reduced 
temperature t (2.50). From (4.26) we extract the (static) critical exponent v. 



^ oc t " oc f 



-1/4 



Z/ 



(4.27) 



Given that the static critical exponent a = |, (4.27) implies that the hyperscaling 
relation is violated 

2- a ^pu, (4.28) 

where p = 3 stands for the number of spatial dimensions of the system. 

Figure 6 shows the inverse of the static susceptibility as a function of q (blue dots) 
right at the critical point k = 2. The solid red line represents the best quadratic (in 
q^) fit to the data 



f -T-critN — 1 
[^fit ) 



-1.57468 ■ 10"® - 1.06109 q^ + 3.84182 q^ + 0{q^) . 



(4.29) 
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Figure 6: (Colour online) The scaling (blue dots) of the inverse of the static suscepti- 
bility xZ=oa ^^ ^^^ critical point, n = 2. The solid red line is a quadratic (in q^) fit to 
the data. 

The red line (4.29) intersects the q^ axis at 

q2 = -1.57468 ■ 10-^ (4.30) 

in excellent agreement with the expected value q^ = (4.10). The data implies 

XZ%,, oc Z"^"* oc q-2 ^^ xZ=o,, ^ q"'^' , (4.31) 

which determines the anomalous critical exponent rj as 

ri = 0. (4.32) 

Figure 7 presents the poles in the dynamical susceptibility at tn = tn* in the vicinity 
of the critical point 

X«i«.,q = 0, (4.33) 

for select values of the momenta q: 

q2 = {10-6,10-5,10-^10"^} ~ {blue , green , orange , black} . (4.34) 

The results of the analysis clearly show that as q — ?■ the values of z% tend to a 
universal profile 

lim i^ = 2.79163 ■ 10"^-0.333392(fi:-2) + 0.0278087(k-2)2 + C((k-2)^) , (4.35) 
q^O q^ 
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Figure 7: (Colour online) Poles of the dynamical susceptibility in the vicinity of the 
critical point, Xreiix),,q = for a set of momenta values q^ = : 10^^ (blue dots) , 10"^ 
(green dots), 10~^ (orange dots) and 10"'^ (black dots). The solid red line is a quadratic 



fit to 2% at q2 



10- 



which is presented by the solid red line on Figure 7^^. Given (4.35) we can determine 
the critical scaling of the relaxation time (see (4.12)) of strongly coupled RN plasma 

(27rTer)-^ =itt), oc q^ ■ (k - kJ oc (27rTc qO^ • (27rT, 0"^ oc (27rTc 0"^ (4-36) 
where we wrote the q dependence as oc q^ and used (4.26). Thus, 

T (X i' (X i^ => z = 4. (4.37) 

In [7] Maeda, Natsuume and Okamura argued that the strongly coupled RN plasma at 
criticality should be identified with the 'model B' according to classification of [8]. As 
such, the dynamical critical exponent z in this model is predicted to be [8] 

z = 4-ri. (4.38) 

Since rj = (4.32) for the strongly coupled RN plasma, we explicitly confirm the 
conclusion of [7]. 

4.3 Universality class of the mass-deformed RN plasma 

In section 2.5 we computed the static critical exponents (a,/3,7,(5) of the conformal 
RN plasma, see (2.53). Further, in previous section we determined the remaining static 



In practice we used the best quadratic fit to to i^ at q^ = 10 
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Figure 8: (Colour online) Poles of the static susceptibility of mass-deformed RN plasma 
in the vicinity of the critical point: Xn=o,s\=<\* ~ 0- '^^^ solid red line is a quadratic fit 
to the data. 

critical exponents (z/, ?]) of this theory, and determined its dynamical critical exponent 
z. 

Notice that the mass deformation of the A/" = 4 SYM plasma modifies its infrared 
(hydrodynamic) properties — it generates a non-vanishing bulk viscosity. Thus, one 
might worry that it is a relevant deformation at criticality and potentially might change 
the universality class of the theory. We show here that this is not the case: the 
universality class of the mass deformed A/" = 4 SYM is the same as that of the conformal 
theory. 

We already argued that some of the static critical exponents of the mass deformed 
theory are unchanged (see section 2.6): 

(«,/3,7,5)=Q,^4'^) ■ ^^-^^^ 

To determine the remaining critical exponents {i>,r],z) we have to generalize the com- 
putation of the dynamical susceptibility as explained in section 4.1 to include the 
fluctuation Zc in addition to {Zh, ^^, Zp}. In analogy to (4.17) we introduce 

Z, = {l-x)-'''' q~^ Zc{x,rv,q), (4.40) 

and solve equations of motion for {zh, -2^, Zp, Zc} with boundary conditions (4.18) sup- 
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Figure 9: (Colour online) The scaling (blue dots) of the inverse of the static suscep- 
tibility of mass deformed RN plasma xZ=oa ^^ ^^^ critical point, k = kc, (4.43). The 
solid red line is a quadratic (in q^) fit to the data. 

plemented with 



lim Zr = finite . 



X— >-l- 



and 



0{x 



3/4N 



as x — 7- Oj 



(4.41) 



As before, up to an overall constant, the dynamical susceptibility is identified with the 
normalizable mode of zj[ near the asymptotic AdS^ boundary, see (4.19). We present 
the results of such analysis for one specific choice of ^, namely. 



M 



10-1x^ = 0.272(1). 



(4.42) 



We use approximate sign since relation (2.59), which was used to obtain (4.42), receives 
corrections of order ^r- 

Figure 8 is equivalent to Figure 5 for the conformal RN plasma. Clearly, the slope of 
the solid red line (representing the quadratic fit to the data) is finite at q^ = 0, which, 
much like in (4.25)-(4.27), implies that the critical exponent i^ = \- Note a technical 
detail: k = k, + 0{M'^/T^); the precise relation is not important for the purpose of 
extracting the critical exponent. From the intersection of the solid red line the k axis 
we find 

kc = 2.0099887(0) . (4.43) 
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Figure 10: (Colour online) Poles of the dynamical susceptibility of mass deformed RN 



plasma in the vicinity of the critical point, Xm= 



ro=ro,,q 



for a set of momenta values 



q^ = : 10 ^ (the solid red line quadratic fit) , 10 ^ (blue dots), and 10" 
The vertical dashed blue line denotes k = kc, see (4.43). 



(green dots) 



Figure 9 is equivalent to Figure 6 for the conformal RN plasma. It represents the 
scaling of the inverse static susceptibility of the mass-deformed RN plasma at criticality. 
As in (4.31), we conclude that 77 = in this case. 

Finally, Figure 10 is equivalent to Figure 7 for the conformal RN plasma. As in 
(4.35)- (4.37) we conclude that z = 4 for the dynamical critical exponent of mass- 
deformed RN plasma. 

5 Conclusions 



In this paper we presented extensive analysis of the critical phenomena in superconfor- 
mal M = 4 SYM gauge theory plasma at finite temperature and a single ^(1) C S't/(4) 
R-symmetry chemical potential. From the explicit analysis of the dynamical suscepti- 
bility near criticality we extracted the static critical exponents (z/, 1]) and identified the 
static universality class of the model^^: 



{a,(3,-f,S,u,ri) 
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This corrects conclusions of [5]. 



Ill 1 

2'2'2'^'4'° 



(5.1) 
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We explicitly computed the dynamical critical exponent of the theory 

z = 4:, (5.2) 

and confirmed the identification of the A/" = 4 SYM plasma dynamical universality 
class with that of 'model B' according to classification of [8], originally made in [7]. 

We demonstrated that although a deformation of the A/" = 4 SYM theory by a 
dimension-3 operator is relevant in the infrared — in particular, it generates the non- 
zero bulk viscosity in the effective hydrodynamic description of the mass-deformed 
plasma — the static and the dynamical universality classes of the theory remains 
unchanged: (5.1) and (5.2). 

We carefully studied the propagation of the sound waves in non-conformal charged 
plasma. We confirmed the computation of the speed of sound from the thermodynamic 
analysis with the direct result extracted from the quasinormal mode analysis of the 
holographic dual. We showed that the bulk viscosity of the mass-deformed A/" = 4 
SYM plasma remains finite at criticality, and satisfies the bulk viscosity bound [13] (at 
least to order O (M'^/T^)). Our computations challenge the Onuki's model [3] for the 
behavior of bulk viscosity near criticality. Indeed, the latter model predicts 



c 



(X \t\-"'+'' (X \t\-^/\ (5.3) 

Onuki 



for the mass-deformed A/" = 4 SYM plasma, in contradiction with the finite result we 
obtained. Thus, when combined with analysis in [2] , it appears that as of now, there is 
no model of transport at continuous phase transitions that is not in confiict with direct 
(first principle) holographic computations. It is important to use gauge theory/string 
theory correspondence to develop a consistent model. 

In order to explain holographic computations of transport at criticality it might be 
necessary to generalize the framework of near-equilibrium relaxation. Specifically, it 
appears necessary to formulate the theory of dynamical critical phenomena in which 
different non-equilibrium correlators relax to equilibrium with different dynamical crit- 
ical exponents. 
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A Coefficients Cij 

X {3k^{x -ly + K+l)- 2H^{5k^{x -lf-K-1)- H^{8k^{x - if - k - 1) 
- 1 - k) . 

(A.l) 

""' 2{x-l){Hl-l){2 + Hl)H,- ^ ■ ^ 

C (2i/g + l)HlK 

'' VT+^{x-l){Hi-l){2 + Hl)- ^ ■ ^ 



2{2Hl + 1))k'HI 
{2 + Hl)m-l){m-l){l + K) 



Ci5 = - .. , ...:;:;;°t;/:/° .. , . • (a.5) 



"-^^ {2 + Hl){Hl-mHl-l){l + >.)- ^ ■ ^ 

_ {2Hl + l){l + K + Hl{{x - ifK^ -2-2k) + Hljl + K))K^H',Pm^ 

i2{i + KnHi-mHi-m2+Hi) ■ ^ ■ ^ 

^- - mm-mlmx-m + .f '''^^' + ^) + ^^ (^^^(^ -'^'-'- '^^ 

+ HI{Ak^{x - 1)2 - 3 - 3/t) + 3 + 3/t) . 

(A.8) 

'- - - ,i2^Him-mx-m^n) ^'''^'^' + ^) + ^^"(^^^(^ - '^' -'^-'^ 

+ 2HI{IIk'^{x - 1)2 - 3 - 3k) + 2Hl{nK^{x - 1)^ + 3fi; + 3) + H^{3k + 3 

(A.9) 
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+ 20fi;2(x-l)2)-3-3/t 



- 3 - Sk) - if^(3 + Sk + 4fi;=^(x - 1)^) - 3 - 3^) . 

(A. 10) 

^^^-" 6(2 + iJ3) • (A-11) 

- 3 - 3k) - iJo^(3 + 3/t + 4/t^(a; - 1)^) - 3 - 3/t) . 



Am 



K 



- H^{32k^{x - 1)2 + 63/t + 63) + 9H^{7 + 7k + Ak^{x - 1)^) + H^^{A5 + 45/t 
+ 32/t2(a; - 1)2) _ H^^AbK + 45 - 2/t2(a; - 1)^)) . 

(A.13) 

/.2gX(l + ^ + gg(^^(^ - 1)2 - 2 - 2k) + Hljl + k)) 

Ak(I-x) 
^31= ^ '-—. (A.15) 

VI + k{H^ - 1) 

6/3/t(l-2/j3)(i_a;) 

C33 - - (,_,)(^/_,).(,^,) (l + - + 4i/3,.(. _ 1). _ 2Hlin + 1 + 2(. - 1)2.2) 
+ 4Hy{x - 1)2 + i7oi2(i ^ ^)) ^ 

(A. 17) 

''" (1 + «:)3/2(i/o' - 1)2(/J3 + 1) • ^^-^^^ 



8k^H^P{x - l)2(-4/73 + 4if6 + 2 + H^^ 
[l + KY/\Hl-l){Hl-l) 



^-35= TTT- ^o/o,.^.Q TTTTTfi T^ • (^A.iyj 



B Sound of A/^ = 4 plasma via the new technique 

The technique for computing the holographic sound wave dispersion relation developed 
in section 3.3 is rather complicated. Thus, we believe that it warrants a simple ex- 
plicit example. This example is being provided by an AdS^ Schwarzschild black hole, 
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holographically dual to strongly coupled A/" = 4 plasma at finite temperature and zero 
chemical potentials. The latter is realized as a special case of the background (2.5) 
with 

A = 0, = 1, x^O. (B.l) 

In this case the only gauge invariant fluctuations are those of zha and zhi. The analog 
of (3.40) is: 

^ " "^^'^ + (1 - x)(x2 - 2a; + 3 - ^2) "^^'^ + x^ - 2a; + 3 - ^2 ^""'^ ' ^ ^^ 

and the analog of (3.42) is 

n - " ^ /9i2 + 3a;^ - 6a; + 1 , 4 

^ " "^^'^ + (l-a;)(a;2-2a; + 3-/3i2) "^^'^ a;^ - 2x + 3 - /3i2 ''''•' ^"^ 

:/h = -^^^ ^., J\ TT-r. il^u + 4/3i2a; - 2/3,2x' - 6/3,2 - 4a;=' + 8/32X 



V3{x - l)(a;2 - 2a; + 3 - /S^y 

v^(x2 - 2x + 3 - /3i2)2 



+ a;^ + 9 + lOa;^ - 12a; - 4/^2 - 4/32X^) 4,o + ^^''^ ^£ ~ ^ll ^T l^\^ ^^ ^h,o • 

(B.3) 



Note that since we expect /3i2 = 1, there are no poles in the connection coefficients in 
(B.2) inside the range of integration x G (0, 1) — the would-be poles are at 

= x^ -2x + 3- (3i2- (B.4) 

Nonetheless, J^h in (B.3) contains factors of (a;^ — 2x + 3 — /3i2)^ in the denominators 
of the connection coefficients. As explained in section 3.3 their origin can be traced 
back to the hydrodynamic expansion (3.25). 

According to discussion around (3.52), solving (B.3) while treating /3i2 as a free 
parameter and imposing the horizon boundary condition as in (3.39) computes Zh,o'- 

l-f3i2 + 2x- x^ 1-Pu .^,. 

ZH,o = ^ — ^ =^ ^o(/3i2) = 7^ — -r- ■ (B.5) 

^ — P12 ^ — P12 

Solving 

^o(/3i2 = /3i2) = =^ /3i*2 = l, (B.6) 

determines the speed of sound c^ = -3^ = |) and following (3.53) computes zh,o- 



ZHfl — ZHfl 



= 2x-x^ . (B.7) 

2o(/9i2)=0 
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Using (B.5) we compute 

= -1. (Bi 



^0 - -7^ 2o{/3i2] 
api2 



2,)(/3i2)=0 



To proceed further, we introduce Zh,i following (3.51). As explained around equation 
(3.55), Zh,i satisfies the following equation 

_ „ 3a;2 - 6a: + 2 , 4 ~ 

° " '"^^ + (l-x)(x2-2x + 2) "^^'1 + x2-2x + 2 ''''' + ^^ 
. 2 , 4 (^-^^ 

^"^ " v^(l-x) '^'° + v^(x^ - 2x + 2) '^'°' 

where we substituted /3i2 = /3i2 = 1- Notice that Jh does not contain factors of 
(x^ — 2x + 2)^ in the denominators of the connection coefficients anymore. Using (B.7) 
we can solve for Zh,i, subject to the horizon boundary condition (3.39): 

Finally, from (3.58) we compute the attenuation coefficient 

/92 = 1, (B.ll) 

which is the expected answer leading to a zero bulk viscosity for the A/" = 4 SYM. 
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